Generation of polarization-entangled photon pairs in a cascade of 
two type-I crystals pumped by femtosecond pulses 
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We report the generation of polarization-entangled photons by femtosecond-pulse-pumped spon- 
taneous parametric down-conversion in a cascade of two type-I crystals. Highly entangled pulsed 
states were obtained by introducing a temporal delay between the two orthogonal polarization 
components of the pump field. They exhibited high-visibility quantum interference and a large 
concurrence value, without the need of post-selection using narrow-bandwidth-spectral filters. The 
results are well explained by the theory which incorporates the space-time dependence of interfering 
two-photon amplitudes if dispersion and birefringence in the crystals are appropriately taken into 
account. Such a pulsed entangled photon well localized in time domain is useful for various quantum 
communication experiments, such as quantum cryptography and quantum teleportation. 
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I. INTRODUCTION 

Photon pairs generated in the process of spontaneous 
parametric down-conversion (SPDC) have been an ef- 
fective and convenient source of two-particle entangled 
states used in tests on the foundations of quantum me- 
chanics as well as application to quantum information 
technologies such as quantum cryptography and quan- 
tum teleportation. In such apphcations, a pulsed source 
of entangled photon pairs is particularly useful because 
the times of emission are known to the users. For exam- 
ple, in the entanglement-based quantum key distribution 
protocol it is convenient to share entangled pho- 

tons between sender and receiver in a train of short op- 
tical pulses synchronized with the common timing clock 
pulses. Such a pulsed source allows the sender and re- 
ceiver to know the arrival times of the photons within 
the duration of the pump pulse and to time stamp each 
received bit readily so that they can sift their raw data 
to generate a shared key from their public discussion. 
A pulsed source is also crucial in a certain class of ex- 
periments that require the use of several photon pairs 
at a time, such as quantum teleportation [Q-^, entan- 
glement distillation [mL and the generation of multipho- 
ton entangled states [p|-p^. In these experiments, time- 
synchronized entangled photon pairs must be available. 
It has also been pointed out that primitive elements of 
a quantum computer can be constructed from the com- 
bined use of several entangled photons and Bell-state 
measurement 111^]. 

A great deal of effort has been devoted to developing 
a pulsed source of entangled photon pairs. Over the past 
decade, femtosecond-pulse-pumped SPDC has been ex- 
tensively studied by several groups. It has been shown 
that the entangled states generated in type-II SPDC 



pumped by ultrashort pulses shows disappointingly low 
visibility of the quantum interference, and narrow- band- 
width filters are required to increase the visibility at the 
expense of the photon flux p^-p^. Theoretical studies 
revealed that the down-converted photons are entangled 
simultaneously in polarization and space-time, or equiv- 
alently, wave number-frequency because of the signifi- 
cant effects of dispersion and birefringence in the crystal 
p3| , p^ 17 1 . As a result of this multiple-entanglement, it 
is possible to distinguish, with some degree of certainty, 
which of the interfering pathways occurs by the measure- 
ment that distinguishes the space-time component of the 
state, for example, the measurement of the arrival time 
of the photons at the detector. This distinguishing space- 
time information is enough to seriously degrade the visi- 
bility of two-photon polarization interference. The pres- 
ence of this distinguishing space-time information stems 
from the fact that the down-converted photons are emit- 
ted spontaneously, but the times of emission are known 
to be within the duration of the pump pulse. The re- 
sultant two-photon state is well localized in space-time, 
and provide undesirable timing capability. The situation 
is serious in the ultrashort-pulse-pumped SPDC, but is 
absent in the ordinary cw-pumped SPDC. Therefore, it 
is vital in the ultrashort-pulse-pumped SPDC to ensure 
that the space-time component of the state contained no 
distinguishing "which-path" information for the photons 
in order to observe the entanglement and quantum in- 
terference in the polarization degree of freedom. Various 
interferometric techniques have been proposed to elimi- 
nate the entanglement in unnecessary degrees of freedom 
and to recover the quantum interference [T7|-|2]]] . 

On the other hand, cw-pumped SPDC exhibits a 
polarization-entanglement with high-visibility quantum 
interference capability. In particular, Kwiat et al. 



1 



realized maximally polarization-entangled photons us- 
ing noncoUinear SPDC with type-II phase matching 
p2[ . More recently, they demonstrated that noncoUinear 
SPDC using two spatially separate type-I nonlinear crys- 
tals pumped by a cw laser exhibits high- visibility quan- 
tum interference |23|. The latter source is particularly 
convenient since the desired polarization-entangled states 
are produced directly out of the nonlinear crystal with 
unprecedented brightness and stability and without the 
need for critical optical alignment. Their result attracted 
great attention for the experimentalists. Recently, Kim 
et al. showed that polarization-entangled photons from 
two spatially separated type-I nonlinear crystals pumped 
by femtosecond laser pulses exhibit high-visibility inter- 
ference [0. In their experiment, coUinear degenerate 
type-I SPDC and a beam splitter were used to create 
a polarization-entangled state. The state prepared af- 
ter the beam splitter was, however, not considered to be 
entangled without amplitude postselection. Only when 
one considers those postselected events (half of the to- 
tal events) in which photons traveled to different out- 
put ports, one can observe quantum interference. Later, 
this problem was solved by using coUinear nondegener- 
ate type-I SPDC and a dichroic beam splitter . It was 
shown that mostly maximally entangled states showing 
high- visibility (92%) quantum interference were obtained 
without amplitude and spectral postselection just before 
detection. However, the available flux of the entangled 
photon pairs was still very limited (<100 s~^). 

In this paper, we apply femtosecond-pulse pumping 
to the second Kwiat scheme with type-I phase-matching 
arrangement. We shall demonstrate that highly entan- 
gled polarization states with sufficiently large flux are 
obtained when an appropriate temporal compensation is 
given in the pump pulses. The results will be shown 
to be consistent with the theory which incorporates the 
space-time dependence of interfering two-photon ampli- 
tudes, if the effects of dispersion and birefringence in the 
two-crystals are appropriately taken into account. 



II. EXPERIMENT 

Our experimental setup is schematically shown in Fig. 
|l|. Two adjacent, thin, type-I crystals (BBO) whose 
optic axes are horizontally (H) and vertically ( V) ori- 
ented, respectively, are pumped by 45°- polarized fem- 
tosecond pulses at 266 nm. The pump pulses were third 
harmonic of a mode-locked output of Ti/Sapphire laser, 
whose approximate average power was 150 mW and rep- 
etition rate was 82 MHz. Due to type-I coupling, H- 
polarized photon pairs at 532 nm are generated by the 
F-polarization component of the pump field in the first 
crystals, and F-polarized photon pairs are generated by 
the ^/-polarization component of the pump field in the 
second crystal. These two possible down-conversion pro- 
cesses equally likely occurs and are coherent with one an- 



other p4| , p5[ . The SPDC was performed under a degen- 
erate and quasi-coUinear condition. The signal and idler 
photons making an angle of 3° with respect to the pump 
laser beam and having the same wavelength around 532 
nm were observed through the identical interference fil- 
ters, centered at 532 nm, placed in front of the detectors. 
The polarization correlations were measured using po- 
larization analyzers, each consisting of a rotatable half- 
wave plate (HWP) and a quarter-wave plate (QWP) (for 
532 nm) followed by a polarizing beam splitter (PBS). 
After passing through adjustable irises, the photons were 
collected using 60-mm-focal-length lenses, and directed 
onto the detectors. The detectors (PMT) were photo- 
multipfiers (HAMAMATSU H7421-40) placed at -1.5 m 
from the crystal, with efficiencies of —40% at 532 nm and 
dark count rates of the order of 80 s'-^. The photodetec- 
tion area was about 5 mm in diameter. The outputs of 
the detectors were recorded using a time interval analyzer 
(YOKOGAWA TA-520), and pulse pairs received within 
a time window of 7 ns were counted as coincident. 

To obtain a truly polarization-entangled state, care 
must be taken to disentangle the polarization degree of 
freedom from any other degrees of freedom, that is, to 
factorize the total state into product of the polarization- 
entangled state and those describing other degrees of free- 
dom. This is equivalent to saying that effective polariza- 
tion entanglement requires the suppression of any dis- 
tinguishing information in the other degrees of freedom 
that can provide potential information about "which- 
polarization" the emitted pair have. In our case, since 
which crystal the origin of each pair is and their polar- 
ization are intrinsically correlated, distinguishing "which- 
crystal" information must be also eliminated. Accord- 
ingly, to make emitted spatial modes for a given pair in- 
distinguishable for the two crystals, we overlap the down- 
conversion light cones spatially by using very thin crys- 
tals each having —130 /im-thickness. In addition, it is 
important to eliminate distinguishing space-time infor- 
mation inherent in the two-photon states produced in 
the ultrashort-pulse-pumped SPDC. Figure ^ schemat- 
ically illustrates what may happen when a 45°- polar- 
ized femtosecond pulse incidents on two cascaded type-I 
BBO crystals, which is estimated roughly from the op- 
tical characteristics of the BBO For simplicity, we 
consider here the degenerate coUinear SPDC. In the first 
crystal, //-polarized down-converted photon wavepack- 
ets are advanced at —74 fs relative to the F-polarized 
pump pulse due to group velocity dispersion, while the H- 
polarized pump pulse is delayed at —61 fs relative to the 
F-polarized pump pulse due to birefringence in the BBO. 
In the second crystal, F-polarized down-converted pho- 
ton wavepackets are advanced at —74 fs relative to the 
//-polarized pump pulse due to group velocity dispersion, 
while the //-polarized down-converted photon wavepack- 
ets are advanced at —242 fs relative to the //-polarized 
pump pulse due to an advance given in the first crystal 
(—135 fs) and dispersion in the second crystal (—107 fs). 
Consequently, after the crystals, space-time components 
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of the two-photon state associated with the polarization 
states \HH) and \VV) are expected to be temporally dis- 
placed by ~168 fs, where \H) {\V)) means a single pho- 
ton linearly polarized along a horizontal (vertical) axis 
and the first (second) letter corresponds to the signal 
(idler). Note that the distance between the two BBO 
crystals is not an important factor in this consideration. 
Since down-converted photon wavepackets are expected 
to have widths comparable to that of pump field (^150 
fs), it is suggested that the space-time components as- 
sociated with \HH) and \VV) do not overlap in space- 
time. As a result, "which-polarization" information may 
be available, in principle, from the arrival time of the 
photons at the detector. 

To eliminate this distinguishing space-time informa- 
tion, a polarization dependent optical delay line for the 
266-nm pump was inserted before the crystals, which is 
denoted as the pre-compensator in Fig. |^. It consists 
of quartz plates, whose optic axes are oriented either 
vertically or horizontally, and a Bereck-type tiltable po- 
larization compensator. The combination of the quartz 
plates of various thicknesses provides a relative delay of 
|T| <350 fs between the H- and F-polarization compo- 
nents of the 266-nm pump field, while the sign is given 
by the orientation of its optic axis. This delay line can 
compensate for relative delay between the state asso- 
ciated with \HH) created in the first crystal relative 
to the state associated with \VV) created in the sec- 
ond crystal to overlap them temporally. The Bereck- 
type polarization compensator is used to adjust the sub- 
wavelength delay. As a result, a two-photon Bell state 
|$+) = {\HH) + \VV))/V2 can be directly created. 

Two kinds of experiments were carried out to con- 
firm whether the target state | <&"'") was successfully pre- 
pared. The first experiment was a quantum state to- 
mography ||2^. The polarization density matrices were 
estimated from 16 kinds of joint projection measurements 
performed on an ensemble of identically prepared photon 
pairs. These joint measurements consist of four kinds 
of projection measurements onto , \V) , \D) , \L)} on 
each member of a photon pair, where \D) = {\H) + 
\V))/V2 and \L) = {\H) + i l"^/V2. We used a max- 
imum likelihood calculation [p8| to estimate the den- 
sity matrix. We also calculated the concurrence of the 
state, which is known to give a good measure for the 
entanglement of a two-qubit system [^. The second 
experiment was a conventional two-photon polarization 
interference experiment. In this experiment, QWPs in 
mode 1 and mode 2 were set to 0°, i.e., their optic axes 
were vertically oriented, while HWP in mode 1 was set 
to ±22.5° and HWP in mode 2 was rotated. In our 
setup, this corresponds to the projection measurements 
onto {\L) or \R)}(E)\e), where \R) = {\H) - i\V))/V2 
and 16*) = cos26' + i sin 261 \V). The visibility of an 
observed interference pattern gives another convenient 
measure for the entanglement. 

At first, we used the interference filters having a full- 



width half maximum (FWHM) bandwidth of 8 nm in 
front of the detectors. Figure |^ illustrates the estimated 
polarization density matrices prepared when (a) no rel- 
ative delay, (b) relative delay of 135 fs, and (c) rela- 
tive delay of 231 fs were introduced between the H- and 
F-polarization components of a 266-nm pump field. In 
these experiments, maximum coincidence counts and typ- 
ical accidental coincidence counts were ~450 cps and <10 
cps, respectively. As shown by these examples, we could 
successfully prepare the states approximately described 
by p{v) = ^{\HH) {HH\ + \VV) {VV\) + v |$+) ($+|. 
These states should exhibit interference in coincidence 
rates for the two detectors in proportion to P{9) — 
j{lzLv sin 46*) . Figure ^ shows the coincidence rates as a 
function of the HWP angle in mode 2 and the evaluated 
visibilities associated with the states in Fig. ^(a)-(c). 
Note that, no background was subtracted in these results, 
so that the evaluated visibilities give the lower limits of 
V. The experiments were repeated for various values of 
delay T introduced in the H- and F-polarization compo- 
nents of the pump field. The evaluated concurrence C 
and visibility v of the polarization interference are plot- 
ted in Fig. I as a function of relative delay T. This 
figure clearly indicates that there is a strong correlation 
between the concurrence and the visibility. This is quite 
natural because C = v should hold for the states p{v) 
given by the above form. In this experiment, the max- 
imum values of the concurrence and visibility were 0.95 
and 0.92±0.02, respectively. The broken line in Fig. ^ 
is a Gaussian curve fitted to the measured data, and the 
FWHM width is 237 fs. 

Next, to investigate the effects of interference filters 
used in front of the detectors, we replaced them with 
those having a FWHM bandwidth of 40 nm. Figure ^ 
shows C and v evaluated by the same procedure, except 
that the backgrounds have been subtracted for both the 
concurrence and the visibility data in this case. The max- 
imum values of C and v were slightly reduced to 0.92 and 
0.90, respectively, while the maximum coincidence counts 
were increased by a factor of « 6. (Unfortunately, the ac- 
cidental coincidence increased more than six times, which 
is probably due to ambient light.) The broken line in Fig. 
^ is a Gaussian curve fitted to the measured data, and 
the FWHM width is «205 fs. The width of the fitted 
curve was slightly reduced in comparison with that for 
the 8-nm filters, which was due to the spectral filtering 
effects, as shown later. 

We also measured the autocorrelation of the pump 
field by measuring the degree of polarization, or 
equivalently, the magnitude of the Stokes parameter 
|s| = \/ s\ + s\ + s\ of the pump just after the pre- 
compensator (and before the crystals), where si, S2, and 
S3 are three independent Stokes parameters. Simple anal- 
ysis shows that |s| is proportional to the temporal overlap 
of the H- and V- polarization components of the pump 
field between which relative delay T was introduced by 
the pre-compensator, and directly gives the autocorrela- 
tion of the pump field. Figure ^ shows the measured |s| 
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as a function of T. We found that |s| can be fitted with 
a Gaussian curve having a FWHM width ^216 fs. As- 
suming the transform hmited Gaussian-shape pulse, the 
FWHM temporal width of the pump field is estimated to 
be ssl53 fs, and that of the pump intensity to be «108 
fs. We see that the width of the autocorrelation of the 
pump field is comparable to that of the measured con- 
currence and visibility data shown in Figs. ^ and ^ This 
result suggests that the width of the measured concur- 
rence and visibility data are determined chiefly by the 
temporal width of the pump field. 



III. THEORETICAL ANALYSIS 

Experimental results are analyzed based on the theory 
given by Grice et al. iQ, Keller et al. and Kim et 
al. ||l^. In the following analysis, we will focus our at- 
tention on the Gaussian curve fitted with the measured 
concurrence and visibility data shown in Figs. ^ and 
^, which will be denoted by V{T). In particular, we 
will analyze how the widths of these curves are deter- 
mined and depend on the experimental parameters. To 
make the following discussion self-contained, it may be 
helpful to begin by showing the minimum of the previ- 
ous theories before discussing our problem. Let us con- 
sider the degenerate, quasi-collinear SPDC in a single 
type-I crystal pumped by a field with the central fre- 
quency Op and the FWHM temporal width t^^"'^^ . 
For simplicity we ignore the transverse components of 
the photon wave vector. Let , Wi] be the frequency 

and kp{ijjp)[ks{ujs), ki{uji)\ the wave number for the pump 
[signal, idler]. The perfect phase-matching condition 
is given by — f2i — $1^/2 for the frequency and 
Ko{^s) + Ko{U,i) — Keiyip) for the wave number, where 
Kj{Q) denotes the wave number for a photon having a 
frequency Q and polarization along the ordinary {j=o) or 
extraordinary (j=e) optic axis. To simplify the notation, 
we introduce detuning for the signal, the idler, and the 
pump from the perfect phase-matching condition by = 
Wj — r2p/2, z/i = Wi — r2p/2, and Vp = Up — flp, respectively. 
To first-order in the interaction, the general expression 
of the two-photon state after the crystal is |14| 
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where a^(z^s) and are the photon creation operators 

for the o-polarized signal and idler modes characterized 
by detuning frequencies Vs and z/^, respectively, which 
are defined after the crystal. The function '^{i/s,i'i) is 
the two-photon amplitude in the frequency domain, and 
is given by Q 



In this expression, a^Vp) is the spectral envelop function 
of the time-dependent classical pump field at the input 



surface of the crystals, 0(i/s, i^i) is the longitudinal phase- 
matching function [|l^, and Cis a constant. Assuming a 
Gaussian shape for the temporal envelop function of the 
pump field, a(fp) is given by 



a(fp) 



(3) 

FWHM bandwidth of the pump field. If z-direction is 
taken to be the pump direction, the longitudinal phase- 
matching function Vi) is given by 



where ap — 



FWHM 

^p 



L/2 



-L/2 



sinc[ ^ L], 



(4) 



where L is the crystal length, A(i^s,i^i) = ks{i^s + ^p/'2) + 
ki{i'i + f^p/2) — kpiup + r^p) is the phase mismatch and 
approximately given by |Tq,n^ 

A(z/3, 1/,) = D+ [us + V,) + -^D" {vs ~ . (5) 

In this expression, D+ and D" are crystal parameters: 
= I/uo(ilp/2) — l/ue(f2p) is the group velocity 
mismatch where Uo{^p/2) {ue{^p)) is the group veloc- 
ity of the o-polarized down-converted photons (the e- 
polarized pump photon) inside the crystal and D" = 



d^R'o/dn^ 



si=o. 



is the group velocity dispersion for the 



down-converted photons. 

Since we are interested in the temporal characteris- 
tics of the down-converted photons, we next consider the 
two-photon state after the crystal in the time domain. 
Formally, it can be written as 



dtsdui>(ts,k)a^tsp{U) |0) 



(6) 



where a^{ts) and fe^(ii) are Fourier transforms of a^(z^s) 
and b^i^i). For example. 



where Lg is the optical path length from the output face 
of the crystal to the observation point. Physically, tg and 
ti are interpreted as the time of arrival of the idler and 
signal at the points that are and Li away from the 
output face of the crystal |Q, and ^{ts,ti) gives infor- 
mation about the space-time motion of the two-photon 
wavepackets. The two-photon amplitude in the time do- 
main is simply given by the Fourier transform of the two- 
photon amplitude in the frequency domain and by |30(| 



(2) ^its,U)^{0\a{t,)b{t,)\-^) 
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dlysdl^,e''^'''''+•'^'^^^ilys,l^^)■ (8) 
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We can further change the variables according to = 
j/s + j/i = Vp, = Vs- Vi, t+ = (ts + ti)/2, and i_ = 
{is — and obtain the equivalent Fourier transform 

relation 



2tt 



X / / dz/+dz/_e-^(''+*++'^-*^V('^+,i^-), (9) 



where 4!{v+,v-) = ^ C^^^*^ i "^ 2" ) V'C^+i*-) = 
+ t+ — Physically, i_ and t+ mean half the 
difference in the arrival time, and the mean arrival time 
of the idler and signal photons, at the points that are Lg 
and Li away from the output face of the crystal, respec- 
tively. The absolute square of the two-photon amplitude 
2 

is proportional to the probability of arrival 

of the idler and signal at the given points for the given 
mean arrival time t+ and half the difference in the arrival 
time Wc show schematically the relevant amplitudes 
in the time domain and the frequency domain in Fig. ^. 

Now, let us proceed to provide a phenomenologi- 
cal model that is consistent with our experiment using 
two crystals. We assume down-conversion processes are 
equally likely to occur in two crystals whose optic axes 
are orthogonally oriented and thy are coherent with one 
another. We represent the overall two-photon state af- 
ter the crystals as a superposition of the states arising 
from the SPDC in each crystal with appropriate phase. 
Then the two-photon state after the crystals in the time 
domain can be written as 



xWffff(i+,t-) \H,t+ + t^)^ \H,t+ - i_>, 
+e'*4'vv{t+,t^) \V,t+ + t^)^ \V,t+ - t_)J (10) 

where the first (second) term represents the two-photon 
state generated in the first (second) crystal, and (p is the 
phase difference between the states generated in the two 
crystals. In Eq. |l^ we incorporated a new dichotomic 
variable in the creation operators and the two-photon 
amplitudes that specifies the polarization {H or V) of 
the photon, and wrote \H,t+ + t^)^ = o,H{t+ + t-) |0) 
etc., which denotes a single-photon state for the sig- 
nal (or idler) mode characterized by a definite polariza- 
tion and time. From the previous discussion, two two- 
photon amplitudes ipHH{t+,t-) and ipvvit+jt-) are ex- 
pected to be temporally displaced by r in the i+-direction 
due to dispersion and birefrengence in the crystals. On 
the other hand, the pre-compensator for the pump be- 
fore the crystals introduces temporal displacement in the 



t+-direction for these amplitudes by —Tp. Therefore, 
we postulate that these two amplitudes can be writ- 
ten using a single common amplitude as ^Phh it+ 7 ) = 
^vv{t++T,t-) = i>it++T,t^), where T = t - Tp. In 
this case, the two-photon state is now entangled both in 
polarization and in space-time, in other words, it is dou- 
bly entangled. The space-time components ij}HH(t+it-) 
and il^vvii+Ti-) provide the distinguishing information 
for the interfering path of the two-photon amplitudes and 
degrade polarization entanglement. The full density ma- 
trix associated with the state given in Eq. |lO| can be 
written as 

Prjkl{ti,t2, h, U) =, (j, t2\g {i, ill*) (*|fc, is), 1^, ti), , 

(11) 

where i,j,k,l — H or V. Since our chief concern is the 
polarization state of the photon pair, we will calculate 
the effective density matrix of the two-photon state af- 
ter the two crystals that are defined in the polarization 
space alone. It can be calculated by partially tracing the 
density matrix in Eq. |l^ over the time- variable, and is 
given by 



(12) 



Pijkl — / / dtidt2Pijkl{tl,t2,ti,t2)- 



Explicit calculation gives the effective density matrices 
of the form 

P=\{\HH) {HH\ + \VV} {VV\ 

+v{T)e-"^ \HH) {VV\ + v*{T)e"^ \VV) {HH\), (13) 

where we used the normalization condition of the two- 
photon amplitude: 



dt+dt^ 



= 1 



(14) 



and notation pijki — p \kl) («, j, k,l = H or V). In Eq. 
[1^ , v{T) is the convolution of the two-photon amplitudes 
in the time domain 



v{T) 



dt+dt-ll!HHit+,t^)tijyy{t+,t^) 

dt+dt^4>{t+ +T,t^)'iJj*{t+,t^). (15) 



It should be noted that v{T) is, in general, a complex 
number. However, in the present experiment, we always 
make i'(T)e~*"^ a real number |f (T)| by using the Bereck- 
type polarization compensator. Then, the density matrix 
given in Eq. ^ is formally the same as what we ob- 
served in the experiment, and Eq. |l^ indicates that their 
concurrence is determined by the convolution of the two- 
photon amplitudes in the time domain. 

Based on the above theory, we have made several nu- 
merical calculations to elucidate the physics behind the 
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measured V{T) curves. We calculated two-photon am- 
plitudes in the frequency and time domains using the 
experimental parameters: Z)+ « —570 fs/mm, D" ft! 855 
fsVmm, L = 0.13 mm, and t^whm ^ ^^53 ^ig- 
ure 1^ shows (a) the (unnormahzcd) two-photon amph- 
tude ^{v-\.,v-) and (b) 'ip{t+,t^). We can see that the 
bandwidth of ^{v+, v-) is much wider in the i/_-dircction 
than in the i^+-direction, while the temporal width of 
ip{t^,t^) is much wider in the t_|_-direction than in the 
i_-direction. This is due to the Fourier-transform rela- 
tions between ■(/)(;/+, ^_) and ?/'(i+,t_). We can also see 
that the temporal width of ip{t^,t^) in the t_|_-direction 
is ft! Tjf^^*^. Using these results, the magnitude of the 
convolution \v{T)\ is calculated and plotted in Fig. |l^ 
together with the fitted curve of autocorrelation of the 
pump field shown in Fig. ^ We can see that \v{T)\ 
nearly agrees with the autocorrelation curve of the pump 
field. This indicates that the temporal shape of V;(i+, t_) 
is approximately Gaussian in the t+-direction and their 
widths are nearly equal to the temporal width of the 
pump field. This is consistent with the fact that, al- 
though there is a small difference in the width, the mea- 
sured V{T) curves nearly agree with the autocorrelation 
curve of the pump field as in Figs. 

Now, let us examine in more detail the origin of the 
different widths of the measured V{T) curves in Figs. || 
and ^ in which the different bandwidth spectral filters 
were used. We will show that this was due to spectral 
filtering effects. So far, the two-photon amplitudes after 
the crystal and before the spectral filter have been dis- 
cussed. They are intrinsic ones in the sense that they 
are determined only by the pump field and the crystal 
parameters and irrelevant to the spectral filters. In con- 
trast, what we actually observed in the experiment was 
the two-photon state after the spectral filters. Hence, it 
is required to relate the state after the filters to the state 
before the filters to find the spectral filtering effects. This 
can be done straightforwardly. Suppose that the linear 
spectral filters used for the signal and idler have complex 
frequency responses expressed by Fs{i^s) and Fi{vi), re- 
spectively. These functions are assumed to be normalized 
such that 



dvi \Fi(yi)Y 



1 



(16) 



holds. Then the two-photon amplitude after the filters 
are simply given by ^'_f(i's, i^i) = '^(vs,Vi)Fs{vs)Fi{vi). 
Since the two-photon amplitude in the time domain is 
related to that in the frequency domain by the Fourier- 
transform relations, the associated two-photon amplitude 
in the time domain is written as 



27r 
1 

(2^ 



dTsdn^{Ts,T^)^(ts - Ts,U - Ti) 



where ^{ts,ti) = Fs{ts)Fi{ti) is a normalized two- 
variable function that incorporates the temporal response 
of the filters, and Fs{ts) and Fi{ti) are the Fourier trans- 
form of Fs{i's) and Fi[ui), for example. 



Fs{ts) 



1 



(18) 



Equation ^ indicates that the two-photon amplitude af- 
ter the filters is given by the convolution of the two- 
photon amplitude before the filters ^'(ts,ti) and the ap- 
paratus function for the filters ^{tg^ti). Because of this 
relation, the input function ^'(is,ii) is smoothed by the 
filter function $(is,ti) to yield broadened output func- 
tion p(ts,ti). It may be worth pointing out that Eq. 
[TtI is formally analogous to the connection of phase-space 
probability distribution with the Wigner functions de- 
scribing system and measuring apparatus (quantum ruler 
or filter) |^,^. In this case, the concurrence and vis- 
ibility are determined by the convolution of the output 
function ipF{t+,t^) = ^ F{(ts + ti)/2, {t^ -ti)/2) in con- 
trast to Eq. |l^, in which the input function 'ijiitj^^t^) 
appears. Accordingly, their widths depend on the tem- 
poral response of the filters, and equivalently, depend on 
the bandwidth of the filters. In general, the smaller the 
filter bandwidth is, the larger the width of the V{T) curve 
is, and vice versa. 

Now, let us discuss quantitatively the difference in the 
width of the measured V{T) curves associated with the 
8-nm and 4G-nm filters. If we assume the frequency re- 
sponse functions Fa{va) and Fi[vi) to be well approxi- 
mated by the Gaussian function with FWHM widths Az^, 
the temporal widths of the associated response functions 
Fs(ts) and Fi(ti) is in proportion to Av~^. In our case, 
Aisnm ~ 104 fs for the 8-nm filter and At4o„m ft 20.8 
fs for the 40-nm filter. On the other hand, from the 
experimental results shown in Figs. || and ^, those con- 
sistent with the present theory are roughly estimated to 

17.2 fs. Here, we 



be At'fi 



8nm 



86.0 fs and At', 



40nm 



used the fact that the width a of the convolution of two 
Gaussian functions each having width ai and a2 is given 
by l/cr^ — l/cr^ + l/cr|. Therefore, even if there remain 
small quantitative discrepancies, it seems reasonable to 
conclude that there is considerable validity in the above 
simplified theory. 

Let us switch our attention to the spectral character- 
istics. For the two-photon state given in Eq. |l|, the 
joint spectral intensity is given hy I {I's, i^i) oc |^'(t^s, t';)!^. 
Then the spectra Is{i^s) and of the signal and the 

idler can be obtained by tracing /(j^s, i^i) over the unob- 
served variable, for example, 



Is{Vs) OC / dVil^L's.l'i). 



(17) 



(19) 



Figure |Tl| shows the calculated individual spectra for the 
down-converted photons just after the crystal. This fig- 
ure indicates that the FWHM bandwidth of the down- 
converted photons is ft 40 nm. Thus, in our experiment 



6 



using the 8-nm filters in front of the detectors, what we 
have observed is only a fraction of the photons having 
the wavelength 532 ± 4 nm that can reach the detectors. 
On the contrary, we can safely say that almost all the 
photons produced in the SPDC have been observed in 
our experiment using the 40-nm filters. In other words, 
the coincidence count rate (« 2700 cps) obtained in the 
experiment using the 40-nm filters had already reached 
the upper limit in our configuration. However, since no 
optimization was made for the transversal spatial mode 
of the down-converted photons in this experiment, we be- 
lieve that further improvement in the coincidence count 
rate may be possible by optimizing the experimental con- 
figuration ||,||. 

While our simplified model seems satisfactory in 
essence, there still remains a problem which needs to be 
solved. Within our model, the autocorrelation curve of 
the pump field should have the smallest width, whereas 
this is not the case. Although we have not arrived at 
a conclusion yet, there might be things neglected in 
this consideration that should actually have been taken 
into consideration. For example, we neglected the ef- 
fect higher than the second-order optical nonlinearity in 
the crystal. There actually might occur higher order ef- 
fects. In particular, it is likely that self-phase modula- 
tion (SPM) due to third-order optical nonlinearity occurs 
when the crystals are pumped by ultrashort optical pulses 
having large peak intensities. When the SPM is taken 
into account, the two- photon state after the crystal is 
chirped, and their coherence length might be shortened, 
and the resultant width of the measured V{T) curves 
might be reduced. To obtain a more satisfactory expla- 
nation, the effect of SPM in the crystals might be incor- 
porated into the theory. 



IV. CONCLUSION 

In conclusion, we have generated pulsed polarization- 
entangled photon pairs by femtosecond-pulse-pumped 
SPDC in a cascade of two type-I crystals. It was found 
that highly entangled photon pairs were successfully ob- 
tained by using giving an appropriate temporal delay 
between the orthogonal polarization components for the 
pump, without the need of spectral post-selection using 
narrow-bandwidth filters. Theoretical analysis showed 
that entanglement depends on the convolution of the 
space-time components of the interfering two-photon am- 
plitudes. The experimental results obtained were well ex- 
plained if the dispersion and birefringence effects in the 
two-crystals are taken into account. Our analysis also re- 
vealed the effect of the spectral filtering on the magnitude 
of the entanglement when the interfering two-photon am- 
plitudes has a different space-time dependence. 

This work was supported by CREST of JST (Japan 
Science and Technology Corporation). 
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FIG. 1. Experimental setup for generating and character- 
izing polarization entangled pulsed photons. The pump pulse 
polarized at 45° irradiates a cascade of two type-I crystals, 
whose optic axes are orthogonally oriented. Pre-compensator 
(quartz plate and tillable wave plate) introduces relative de- 
lays between the two-photon amplitudes created in the first 
and second crystal. CC is the coincidence circuit. 
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FIG. 2. Schematic of the temporal development of the 
down-converted photons and the pump. See text for details. 
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(a) r=Ofs C=0.404 (b)r=135fs C=0.956 (c) 7"=231 fs C=0.651 




FIG. 3. Estimated polarization density matrices for (a) no relative delay, (b) relative delay of 135 fs, and (c) relative delay 
of 231 fs between the H- and V-polarization components of the pump field. Only the real parts of the matrices are shown. The 
contribution of the imaginary parts is negligible. The dark-colored bars show diagonal elements, while the light-colored bars 
show off-diagonal elements. Estimated values of concurrence are shown as C. 
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FIG. 4. Polarization correlation experiment. Coincidence counts in 10 s is plotted against the orientation angle of the 
half- wave plate in mode 2, while the half- wave plate in mode 1 was fixed to ±25.5°, for the states shown in (a)-(c) of Fig. ^. 
Estimated visibilities are also shown in the figure. 
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FIG. 5. Evaluated concurrence C and visibility v when 
8-nm interference filters were used. These values are plotted 
against the relative delay introduced between two orthogonal 
components of the pump field. The broken line is Gaussian 
fit to the data with 237-fs widths (FWHM). 
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FIG. 6. Evaluated concurrence C and visibility v when 
40-nm interference filters were used. The broken line is Gaus- 
sian fit to the data with 205-fs widths (FWHM). 
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FIG. 7. Autocorrelation of the pump field evalu- 
ated by measuring the degree of polarization after the 

prc-compcnsator. They are plotted against the relative delay 
introduced between two orthogonal components of the pump 
field. The broken line is Gaussian fit to the data with 216-fs 
widths (FWHM). 




FIG. 8. Illustration of two-photon amplitudes in the time 
and frequency domains generated by the femtosecond-pulse 
pump. For simplicity the case where Ls = Li = L is dipicted. 
They are connected by a Fourier transform relation. 
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FIG. 9. Calculated (unnormalized) two-photon amplitudes (a) in the frequency domain and (b) in the time domain. Only 
the real part of the amplitude is shown for the time domain. The contribution of the imaginary part is relatively lower. 
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FIG. 10. Calculated magnitude of the convolution |«(r)|. 
The broken line shows the autocorrelation of the pump field. 
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FIG. 11. Calculated spectrum for the down-converted pho- 
tons just after the crystal. Slight assymetry in the spectrum 
with respect to the phase-matching wavelength (=532 nm) is 
due to the combined effect of the group velocity mismatch 
and the group velocity dispersion in the nonlinear crystal. 
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